Electronic Raman scattering in a multiband model for cuprate superconductors 
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Charge-charge, current-current and Raman correlation functions are derived in a consistent way 
using the unified response theory. The theory is based on the improved description of the conduction 
electron coupling to the external electromagnetic fields, distinguishing further the direct and indirect 
(assisted) scattering on the quasi-static disorder. The two scattering channels are distinguished in 
terms of the energy and momentum conservation laws. The theory is illustrated on the Emery three- 
band model for the normal state of the underdoped high-T c cuprates which includes the incoherent 
electron scattering on the disorder associated with the quasi-static fluctuations around the static 
antiferromagnetic (AF) ordering, ft is shown, for the first time consistently, that the incoherent 
indirect processes dominate the low-frequency part of the Raman spectra, while the long-range 
screening which is dynamic removes the long-range forces in the A\ g channel. In the mid-infrared 
frequency range the coherent AF processes are dominant. In contrast to the nonresonant B\ g 
response, which is large by itself, the resonant interband transitions enhance both the Ai g and 
Big Raman spectra to comparable values, in good agreement with experimental observation. It 
is further argued that the AF correlations give rise to the mid-infrared peak in the Bi g Raman 
spectrum, accompanied by a similar peak in the optical conductivity. The doping behavior of these 
peaks is shown to be correlated with the linear doping dependence of the Hall number, as observed 
in all underdoped high-T c compounds. 

PACS numbers: 78.30.-j, 74.72. Dn, 74.25.Gz 

Keywords: high-T c superconductors, electronic band structure, long-range Coulomb screening, electronic 
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I. INTRODUCTION 



Multiband models often present several energy scales 
of the same order of magnitude, related to various anti- 
crossings of the bands. One such interesting example is 
the Emery model for the high-T c cuprates. The effective 
band structure of this model exhibits hybridization gaps 
related to the anticrossings of three bands associated with 
the Cu02 unit cell, as well as the dimerization pseudo- 
gaps related to the antiferromagnetic (AF) fluctuations, 
all of the order of 0.1 eV. The obvious prerequisite for the 
understanding of the high-T c superconductivity, which in 
turn is associated with energies of the order of 0.01 eV, 
is the correct identification of the origin of the 0.1 eV 
energy scales. In the attempt to distinguish among the 
0.1 eV energy scales, one is left only with the difference 
in the associated behaviors in the momentum space, i.e. 
with the corresponding coherence factors, to use the band 
language. As is well known, the coherence factors reflect 
the crystal symmetry and experimental probes sensitive 
to the associated selection rules, such as infrared conduc- 
tivity and Raman scattering [E S S 0, S, @, 0, 1, i] , are 
well suited for the study of the coherence factors. The 
motivation of the present paper is to discuss theoreti- 
cally the existing Raman data from such a point of view. 
This is accompanied here by the solution of several long- 
standing problems which concern the electronic Raman 
scattering in general. 

More specifically, the experimental Raman investiga- 
tions of the effects of superconductivity on the Drudc 
part of the B 2g spectra of YBa 2 Cu307_ a: 0] and 
Bi2Sr2CaiCu208+a; [H| confirmed the conclusions of 



other experiments [ill Il2j that the superconducting 
gap/pseudogap is of the order of 25 meV, with a predom- 
inant d x 2_ y 2 symmetry. In addition, the B\ g spectra in 
underdoped I^-zSr^CuC^ [l3[ and Bi2Sr2CaiCu208+a; 
Q compounds show at temperatures up to room tem- 
perature a strong two-magnon peak at 0.1-0.3 eV and 
a secondary structure at a frequency about three times 
lower. Both scales exhibit the same doping behavior. 
The smaller scale is therefore usually associated with the 
single-paramagnon AF pseudogap [f| H, [TH, [r| ■ Similar 
scales appear in other experiments, in particular in mea- 
surements of the specific heat 12j. Equally important 



are the overdoped cuprates where 0.1 eV energy scales 
are observed in featureless mid-infrared spectra in op- 
tical conductivity and Raman experiments [H, 0> S @1- 
The latter are usually associated with the strong quasi- 
particle damping effects, that is, with the scattering from 
the uncorrelated spin disorder, rather than with the AF 
paramagnons and the concomitant disorder. 

The small energy scales of the order of 0.1 eV and 
less occur in the Emery three-band model for the high- 
Tc cuprates in the limit of large repulsive interaction on 
the Cu site [ijj]. This interaction is renormalized out by 
introducing the auxiliary bosons [l||, which forbid the 
double occupancy of the Cu site, i.e. by introducing the 
Mott charge correlations. The result for finite doping 
is the effective band structure with bands broadened by 
the scattering of fermions on bosons. The single-particle 
dispersions obtained on the hole-doped side within the 
paramagnetic non-crossing-approximation (NCA) or 
dynamical mean-field theory [171 ] approaches are simi- 
lar to those found by the simple mean-field slave-boson 
(MFSB) theory [l5|, [l8j], when the latter is supplemented 
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by harmonic boson fluctuations around the mean-field 
saddle-point. The band dispersions introduce the non- 
magnetic energy scales of the order of 0.1 eV and less, 
in particular through the splitting between the resonant 
band and the main band. The band broadening r(k, ui) 
of the non-Fcrmi-liquid type is related to the inelas- 
tic scattering on anharmonic (slave) bosons, which de- 
scribe the Cu-0 charge fluctuations irrespective of the 
spin. r(k, w) is itself characterized by the energy scales 
of 0.1 eV. The Raman background corresponding to the 
charge fluctuations was evaluated within the NCA [1(1 . 
It reflects the same non-magnetic 0.1 eV energy scales, 
in particular through the processes of charge excitations 
from the main oxygen band to the resonant band. The 
agreement between the calculated single-particle [l8[ and 
electron-hole Raman II 611 pr operties and the correspond- 
ing ARPES [H, EH En^l and Raman 0, [| measure- 
ments on La2- a: Sr a ;Cu04 family of materials is remark- 
able. 

In this kind of approach the magnetic effects manifest 
as perturbations in terms of AF paramagnons [23f ]. The 
associated pseudogap energy Aaf is well below 0.1 eV. 
Until now, the bosonic effects of paramagnons were es- 
timated only by omitting the band broadening due to 
bosonic charge fluctuations. This amounts to the use of 
the MFSB theory, supplemented by the coupling of the 
Fermi liquid to the paramagnons [23j ]. Such an approxi- 
mation conserves the 0.1 eV energy scales in the band dis- 
persion and allows for the (in)elastic scattering on para- 
magnons. The corresponding inelastic processes turn out 
to be more important [23[ on the hole-doped side than 
on the electron-doped side of the "non-magnetic normal 
state" extrapolated close below the superconducting T c . 
The whole hierarchy of energy scales, and especially the 
assertion that the relevant non-magnetic energy scales 
are larger than Aaf, which itself is larger than T c , is ob- 
viously of essential importance for the understanding of 
high-T c superconductivity. 

In order to investigate carefully the energy scale hi- 
erarchy, this paper is focussed on the effect of the AF 
paramagnons on the Raman response, introducing fur- 
ther simplifications which nevertheless conserve the main 
non-magnetic and magnetic scales at and below 0.1 eV. 
The nonmagnetic scales below 0.1 eV are retained in the 
fermion dispersion. The AF correlations are described 
by the AF gap Aaf instead of the pseudogap [24| and 
by the bosonic fluctuations (magnons) around the AF 
state. Both steps are usually considered as legitimate 
for temperatures below Aaf [HI]. In this way, the in- 
elastic scattering on magnons is neglected (in addition to 
that on charge fluctuations). This omits in particular the 
antiadiabatic magnon effects on the single-particle spec- 
trum of holes [23| at energies very much below Aaf ■ The 
whole approach reduces in this way to the MFSB three- 
band theory with the AF dimerization which includes 
(only) the elastic scattering on the (intrinsic AF and ex- 
trinsic) disorder. Even with such drastic simplifications 
the problem is a serious one. 



This article investigates in detail the Raman spectra 
of the underdoped cuprates and distinguishes among the 
coherence factors associated in the reciprocal space with 
the non-magnetic and magnetic scales which appear in 
the problem. The usual Raman analysis of the high-T c 
cuprates starts from the simple Abrikosov and Genkin 
approach 26]. The latter treats the bi- linear Raman 
excitations as non-resonant and calculates the Raman 
intraband contributions starting from the free electron 
limit 0, SI SI SI M, E2, El SI 0, S|. This is 
replaced here by the description of the electron-photon 
coupling effects which is more appropriate for the analy- 
sis of the relevant coherence factors for a nearly half-filled 
tight-binding band. In such a discussion it is obviously 
important to account also for the decoherence effects, 
associated at least with the elastic scattering of charge 
carriers on the quasi-static disorder. 

Associated is the problem of the screening of the long- 
range Coulomb forces in the presence of the disorder 

Mil Si, M EH, El, SI Si EE SE S3, Si, [p- This 

problem is usually treated in the Raman (and infrared) 
analysis by the field-theory approximation (FTA). In this 
approach the long-range forces are screened off by the co- 
herent long-range screening and the elastic scattering on 
the disorder is taken to break the translational symme- 
try, i.e. the momentum conservation laws. The two steps 
may thus seem to be either contradictory or to amount 
to double counting. By distinguishing the (direct) pro- 
cesses with the quasi- particle momentum conservation, 
from the (indirect) processes, which do not conserve the 
momentum, we show therefore that the two steps in ques- 
tion can be reconciled. The momentum conservation pro- 
cesses are subject to the coherent long-range screening, 
while the other processes do not imply long-range forces 
at all. 

Being interested here primarily in the interband scales 
we extend the above single-band considerations to the 
multiband case. The role of interband transitions is 
twofold here. First, the quasi-particles can be excited 
resonantly from the conduction band to the other bands. 
Second, the excited quasi-particles relax back into the 
conduction band, assisted by the elastic scattering on 
the disorder. The former effect is treated by replacing the 
usual static-Raman- vertex approximation (SRVA) by the 
elastic- Raman- vertex approximation (ERVA). This rep- 
resents a natural extension of the recent multiband op- 
tical conductivity analysis [13, [4l[ to the Raman case. 
Such an approach gives access to the most important 
non-magnetic single-particle scales of the Emery model. 
On the other hand, it is shown that the additional elastic 
scattering on the disorder, associated with the interband 
transitions, can be included into the (indirect) processes, 
which do not conserve the quasi-particle momentum. 

The result of these steps is the theory of the electronic 
Raman scattering in multiband models, the Emery model 
for the high-T c cuprates in particular, which can be com- 
pared to the experimental findings. As the analogous 
theory applies also to the conductivity, this approach, 
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as a whole, establishes the relation among a number of 
measurable quantities including the DC conductivity and 
the Hall number, all sensitive to the anomalous features 
in the quasi-particle spectrum close to the Fermi level, 
such as hybridization/dimerization (pseudo)gaps and the 
van Hove singularities. It appears that the AF dimer- 
ization gap produces the intensity maximum in the B\ g 
Raman channel as well as in the optical conductivity, 
while the low-lying E>2 g spectrum remains unaffected. 
In addition, the number of the van Hove singularities is 
doubled, which restores approximately the local electron- 
hole symmetry in the conduction band. This agrees fully 
with the measured doping dependence of the Hall num- 
ber in the underdoped electron- and hole-doped regimes 
[3, The small 0.1 eV energy scale observed in all 

these experiments in the underdoped cuprates is thus as- 
sociated here with the AF dimerization rather than with 
the non-magnetic scales of the same order of magnitude. 
Such interpretation requires however further confirma- 
tion through the theory beyond the MFSB level. 

The paper is organized as follows. In Sec. II the re- 
sponse of the electronic system to external transverse 
vector fields is formulated for a multiband model and 
applied to the Emery three-band model where the lo- 
cal field corrections are absent. The contributions of the 
direct and indirect electron-hole pair excitations to the 
Raman correlation functions are determined, including 
the screening by the multiband RPA (random phase ap- 
proximation) dielectric function. The structure of the 
low-frequency (Drude) contribution to the Raman corre- 
lation functions is given in Sees. Ill and IV. The relation 
between the ERVA and SRVA is discussed in Sec. V. The 
predictions of the model with AF correlations for the Hall 
number, the optical conductivity, and the correspond- 
ing contributions to the B\ g and i?2g Raman spectra are 
given in Sec. VI, and compared to the experimental data. 
Sec. VII contains the concluding remarks. 



II. MULTIBAND MODEL HAMILTONIAN 
A. Emery three-band model 

We consider the conduction electrons described by 
the reduced version of the quasi-two-dimensional Emery 
three-band model [ijj], in which the second- neighbor 
bond energy t pp is set to zero, and the short-range inter- 
actions V p d and V pp are approximately included in the 
copper and oxygen single-particle energies. The Hamil- 
tonian is 

H = H Q + H[ + H' 2 + H cxt . (1) 

Ho is the effective single-particle term. The electron 
quasi-elastic scattering on the disorder is described by 
H[. H' 2 = H c + Haf represents the two-particle inter- 
actions, including both the long-range Coulomb forces 
(H c ) and the residual interactions responsible for the AF 



correlations (Haf)- _ff oxt describes the coupling of the 
conduction electrons to the external fields. 

Using the slave-boson approach to treat the limit of 
large Hubbard interaction on the copper site Ud, the 
effective MFSB single-particle Hamiltonian [l5[ can be 
written in the representation of the non-diagonal trans- 
lationally invariant states as 

Ho = E [<(kCC + H.c.], (2) 

U'ka 

with the orbital index l,V — d, p Xl p y . Here the di- 
agonal and off-diagonal matrix elements have the well- 
known form: Hf} (k) = Ei — 2tj^cosk z a 3 , i?g P "(k) = 
2it^sin|k-a a , with a = x,y, and H^ Py {k) = 
— 4i pp sin 5k • ai sin ^k • &2 (ai, a2, and a3 are the primi- 
tive vectors of the tetragonal lattice in question) . Ei are 
the renormalized site energies, t p ^ is the renormalized 
first-neighbor bond-energy, t pp is the second-neighbor 
bond-energy, and t± is the interplane bond-energy. Using 
the transformations 

L 

Ho is diagonalized in terms of three bands 

Ho = E^( k )4 CT ^ 5 (4) 

ilea 

with the band indices L — c for the nearly half filled 
(conduction) bonding band and L — N, P for the non- 
bonding and antibonding bands (which are empty in the 
hole picture used here). For t pp = the structure of 
E L (k) and U w (l,L) is well known [HEH]. 

The effects of the AF correlations on the Raman spec- 
tral functions is approximated here by replacing the cou- 
pling of the conduction band electrons to the AF fluctua- 
tions by their coupling to the Qaf mode, which is taken 
as frozen in. The effect of bosons with the wave vectors 
close to Qaf on the quasi-particle dispersion is thus ne- 
glected, i.e. the pseudogap is replaced by the gap A(k) 
involved in Haf US [Hf, 

Haf = E [A(k)4 ff c k ±Q AF , + H.c] . (5) 

kcr 

On the other hand, the life-time effects associated with 
slow AF fluctuations can be included in the H[ quasi- 
elastic scattering on the disorder 0, 5i| , 

H[ = E ^"(k-k'^LW (6) 

Lkk'tr 

This implies the adiabatic approximation in the quasi- 
particle scattering on bosons, i.e. the boson frequency 
lower than the temperature of interest [23j]. As already 
pointed out in Introduction, the corresponding correc- 
tions are not expected to affect much the conclusions 
which concern the 0.1 eV scale in the underdoped com- 
pounds, below the two-magnon resonance 0, Q. This is 
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the range to which we restrict ourselves here, while dis- 
cussing some basic questions, which concern the Raman 
scattering itself. 

Finally, the long-range forces are given by 



Hc = El^(-q)?(q), 



q/0 



vq 



(7) 



with q(q) being the charge density operator, 

= EE^ iL '(k,k+q)4x +qc n (8) 



LL' ko 



and the q LL (k, k + q) are the related dimensionlcss 
intra- and interband charge vertices [see Appendix C and 
Eq. ©]. 



B. Electromagnetic coupling 

The coupling of the conduction electrons to the elec- 
tromagnetic fields polarized in the a and/or /3 direction 
follows from the minimal gauge-invariant substitution 

(Slum 



ff ext = F ext + ff ext = _I V] A Q (q)J Q (-q) 



2mc 2 



A a (q-q')^(q')7^(-q;2). (9) 



qq'a/3 



Here 



U<i) = EE J "'( k )44 +qffl 

LL' kcr 

7 Q/3 (q;2) = ^^ 7 ^'(k; 2)4X +qtT , (10) 

LL' kcr 

are, respectively, the current density and bare Raman 
density operators [ID, H(| ■ The explicit form of the cur- 
rent vertices, J^ L (k), and the bare Raman vertices, 
7^3 (k; 2) for the i pp = Emery three-band model are 
given in Appendix A. 

The coupling (|9|) can be completed with the coupling 
to the external scalar fields F GXt (q), 



Hf* = ^F ext (q)g(-q), 



(11) 



used in the longitudinal response theory (see Appendix 
C). It is important to notice that, due to the absence of 
the local field corrections [H,|4^| in the Emery model, the 
long-wavelength charge vertices (q = J2 a 1ae a is small) 
satisfy the general relation [H, IH3| 



eq L ' L {k + q,k) « e8 L ,L> 

+ (1-^lOE 



(12) 



hq a jL' L (k) 



III. RAMAN CORRELATION FUNCTIONS IN 
PURE SYSTEMS 



In the mean-field slave-boson theory [l5j used here, the 
physical Raman correlation functions are proportional to 
the corresponding correlation functions of the auxiliary 
fermions described by the band structure associated with 
Eqs. dU) and It goes without saying that the same 
conclusions hold for the physical fermions with the negli- 
gible local interactions Ud- The simplest operative way to 
determine the Raman correlation functions of this three- 
band auxiliary fermion model is to consider the Gold- 
stone theorem for the thermodynamic potential in the 
Matsubara representation with H' = H cxt + H' 2 + H[ rep- 
resenting the perturbation, and collect all fourth-order 
contributions in the vector fields ^4 a (q") and Ap{q[). It 
is convenient to divide this procedure into four steps. 
First, the H' = i/ cxt case provides the definition of the 
Raman vertex functions in the multiband model under 
consideration, with particular care devoted to the res- 
onant enhancement of the Raman scattering processes. 
Second, for H' = H ext + H c , we shall define the direct 
contributions to the Raman correlation functions and re- 
consider the role of the long-range screening in the pure 
multiband models. Third, by considering the perturba- 
tion H' = H ext + H c + H[, we shall introduce the distinc- 
tion between the direct and indirect (disorder-assisted) 
electron-hole excitations and discuss which of these pro- 
cesses dominate the Raman spectra measured in the high- 
T c cuprates. Finally, by including Haf, we shall study 
the influence of the low-frequency excitations across the 
AF (pseudo)gap on both the Drude part and the related 
low-lying interband part of the Raman spectrum. 



A. Raman vertex functions in pure systems 

In the absence of the disorder and AF scattering pro- 
cesses, the direct summation of the fourth-order diagrams 
in the vector fields A a (q") and ^(q') leads to Fig. 1(a), 
representing the Raman correlation function in the ideal 
lattice, approximately given by its intraband contribu- 
tion. Namely, in the high-T c cuprates, the interband exci- 
tation energies are of the order of typical optical energies, 
1.75-2.75 eV, which is far above the largest Raman shift 
(defined below) measured in experiments (hui < 1 eV). 
Consequently, the interband contributions to the Raman 
correlation functions can safely be neglected in the ideal 
lattice. As will be seen below, the AF correlations intro- 
duce the possibility of the low-lying "interband" excita- 
tions requiring the generalization (Sec. Ill B) of Fig. 1(a). 

Thus, in a pure system (denoted by p) we have 



(13) 



with the longitudinal current vertices L (k) identical 
to the transverse current vertices given by Eqs. (JTT 



kk', 



x -X>£ c (k, k+, k' + , k', uA 7 ^(k', u „ 



5 



(a) 




(b) 



FIG. 1: (a) The purely electronic intraband Raman correla- 
tion functions in a pure system, (b) The Raman vertex (full 
rectangle) shown in terms of the bare Raman vertex (full cir- 
cle) and the interband current vertices (open circles). 



where T>p C (k, k + , k^, k', ui) is the intraband electron-hole 
propagator in the ideal lattice, defined by 



-V 



p LL '(k,k+,k;,k» (14) 
/ L (k)-/ L ,(k + q) 



Jk,k' 



hiu + E L (k) -E L >(k + q) 



vq 



for the band indices L = L' = c. /x(k) = /(£x(k)) is 
the Fermi-Dirac distribution function. Furthermore, the 
7^g(k, u>i, u) s ) are the related intraband Raman vertices 



7^0'^-^) = — 2 

J^(k)J^(k) 



J Q ic (k)J^(k) 



huj. 



fuv s + ElcQz) + if] 



- If] 



(15) 



and k + = k + q. Here u>i,q",a and w s ,q',/3 are the 
frequencies, wave vectors and polarization indices of the 
incoming and scattered photons, respectively. u> = u>i — 
u! s is the Raman shift, q = q" q', and Ell'O^) — 
El (Is.) — El'Qs)- Eq. (fT5")) is gauge invariant in the limit 
77 — ► 0. As mentioned at the beginning of this section, 
both the scattering processes on the disorder and the AF 
correlations are absent in 7^3 (k, u)i, lo s ). 

The diagrammatic representation of the Raman ver- 
tices is shown in Fig. 1(b). The first term on the right- 
hand side is the quadratic coupling term, while the latter 
two represent the bi-linear contributions. The resonant 
nature of the Raman scattering processes refers to the 
bi-linear terms. The resonant effects are large in the 
high-T c cuprates because, as mentioned above, the in- 
terband excitation energies ElcQ*) m Eq. (fT5"]) are of the 
order of typical optical energies. In addition to the res- 
onant condition, £x c (k) w hiOi and/or £x c (k) « huj s , 
the efficiency of the resonant enhancement of the Ra- 
man scattering processes depends also on the relaxation 
processes in the intermediate interband photon absorp- 
tions/emissions that are omitted here. Although, in prin- 
ciple, these relaxation processes have to be treated on 
an equal footing with the relaxation processes in the 



electron-hole propagators V LL (k, k + , k+, k', u>), we shall 
use below an approximate treatment, by including the 
former phenomenologically (see Sec. Ill A. 2) and the lat- 
ter by using the direct summation method (Sec. IV). 



1. Effective mass theorem 

Let us consider the LOi = lo s = 0, 7/ — > limit of 
Eq. (fT5")) . The result is the static Raman vertex of the 
form 



l c a %(k) = ^Jk-,2) + - 



m ^ 2j£ c (k)Jf (k) 



2 Z-^i 



(16) 



Here the symmetry relation J^ c (k) = J^ L (k) has been 
used. This expression can be combined with the relation 



7S( k ) 



m d 2 E c (k) 



to obtain the "effective mass" theorem 
m d 2 E c (k) 



h 2 dk a dkp 



7^(k;2) 



^2 



2J^(k)JI L (k) 
E Lc {k) 



(17) 



Eq. (fT7|) [and Eq. (fl~9|) ] holds even when its left-hand 
side is dependent on k, i.e. beyond the effective mass 
approximation in the vicinity of the Fermi level. The 
result is appropriate for any multiband model with the 
hole-like (— sign, the case considered here) or electron- 
like (+ sign) dispersion of the conduction electrons. 

Eq. (fTT)) turns out to be important for both 
the conductivity-sum-rule analyses and the transport- 
coefficient studies, in particular when the AF term ([5]) is 
included. Actually, Eq. (fTT|) represents a partial conduc- 
tivity sum rule for three bands [40], which holds when 
the photon frequencies are small with respect to the 
transition frequencies into all other bands. When the 
high-frequency transitions are included in the present ap- 
proach the "effective mass" is replaced by the free car- 
rier mass, i.e. the present tight-binding (Wannier) ap- 
proach [43| satisfies the general sum rule established by 
Abrikosov and Genkin [2g, [27| • 

The theorem states that the zero-frequency electron- 
hole pairs (corresponding to the formal limit LUi,Lu s — > 0) 
can be excited by the electromagnetic fields through 
the bare quadratic electron-photon coupling and/or 
through the bi-linear term in which the first-order (high- 
frequency) interband excitations appear as virtual inter- 
mediate states. 



2. Elastic- Raman-vertex approximation 

Since the Raman shift uj = — u s is small in compar- 
ison with the typical values of u>i or u> s , it is reasonable, 
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in the numerical calculation in Sec. V, to use the elastic- 
Raman- vertex approximation 

7^(k,cJi,w s ) w 7^(k,w i ,a; i ) = 7^(k,u>i), (18) 

in which the zero-frequency processes (oji , oj s ~ 0) are ap- 
proximately separated from the higher-frequency absorp- 
tion/emission processes. The phenomenological treat- 
ment of the interband relaxation processes in the res- 
onant channel then gives rise to the general gauge- 
invariant expression which reduces to Eq. (|15|) in the limit 

r intor /wi -> o 



(^) 2 J^(k)jf (k) 

££ c (k) 



x ■ 



2£ ic (k) 



(ft* + ihT inter ) 2 - #£ c (k) 



(19) 



[again J^ c (k) = J£ L (k) is used]. 

It is useful now to incorporate the symmetry properties 
of the Emery three-band model into Eqs. (fT5|) and (JTHJ). 
First, we remember that the analysis of the electronic 
Raman spectra of the high-T c cuprates is usually focussed 
on the in-plane polarization of the electromagnetic fields 
(a, f3 = x,y). It is thus convenient to arrange the Raman 
vertices according to the irreducible representations of 
the D 4 k point group [U H3, [ElJ . The resulting Raman 
vertices are of the form 7^ c (k, oji), with the label v — A\ g , 
Big, and B 2g representing the Ai g , B lg , and B 2g Raman 
channels, respectively. The symmetrized vertices are 



7i c l9 (k^i) 
7B 2s (k,^) 



7^(k,wi) -7^(k,cJi), 
T£(k,wi). 



(20) 



It should be noticed here that the Raman correlation 
functions of the tetragonal high-T c cuprates are diago- 
nal in this representation. The orthorhombic distortion 
of the Cu02 plane, which occurs in some compounds 
(YBa2Cu307_a;, for example), mixes these three chan- 
nels. However, as previously estimated [43j], the mixing 
is typically of the order of 1/10 and is neglected in the 
present analysis. 



B. Long-range screening in pure systems 

The effects of the long-range Coulomb forces on the 
Raman correlation functions are given in the usual way 
HHSHSSIil. In absence of the incoherent 
scattering processes those functions are described by the 
diagrams in Fig. 2(b). The screened correlation function 
Xv,v{<\,u,Ui) is given by 



(21) 



47re 2 

-x*,i(q 5 oj, Ui) ^* > xi,f(q> w > 




(a) 



k+q k+q k+q k+q 

'v k \ \ k \ X k k X 

FIG. 2: (a) The Raman correlation functions in a general 
case with the long-range forces and the quasi-elastic scattering 
processes taken into account. The full rectangle is the Raman 
vertex of Fig. 1(b). The shaded box includes the electron-hole 
self-energy contributions associated with both the long-range 
forces and the scattering processes on the disorder, (b) The 
long-range screening of the Raman correlation functions in 
the case where the scattering processes on the disorder are 
absent. The open circles represent the charge vertices and 
the dashed line is the long-range force Awe 2 /q 2 . 



The coupling function x„,i(q, oj, oj{) is defined by 
Eq. dm), with 7^(k,w-) 7 ^(k / ,a; i ) replaced by 
7^ c (k, uji)q cc (k' + q, k'). The dielectric function in 
Eq. (|2"Tj) has the form 



e(q,oj) = £oo(q,w) 



4-7re" 



Xi,i(q,w), (22) 



with e 2 xi.i(q, oj) representing the charge-charge correla- 
tion function given by 



3 2 XM ( q ,c) = i]T5>V L '(k,k + q) 



(23) 



LL' kk'tr 



xq LL (k' + q, k') -V LL (k, k + , k' + , k' , oj) . 

Here V LL (k, k + , k' + , k', u>) is the electron- hole propaga- 
tor defined in Appendix C. 

For the B\ g and B 2g Raman channels, the coupling 
functions Xv,i(q, oj , u>i) vanish for symmetry reasons, and 
the long-range forces do not affect the Raman spectra in 
the Bi g and B 2g channels. Furthermore, it is useful to 
separate the constant term in the A\ g Raman vertex from 
the dispersive term (2(| [2?j|, 7^ c lg (k, ojj) — 7% (Wj) + 
7A lg (k,^i)i in tne wa Y tnat XAi g ,i(q ; oj, oji) = [no- 
tice that 7^ c (k,Lj i ) = 7° c (k,Wi) for v = B lg ,B 2g }. In 
this way Xf,i(qj w j w i) = f° r a ^ three Raman chan- 
nels. [The hat in Xv.i{<\,oJ,oji) indicates that only the 
dispersive part of the vertex 7^ c (k, u>i), 7^ c (k, w,), is in- 
cluded in Xv,i(<l,0J,0Ji).] Consequently, the dispersive 
terms 7^ c (k, Wj) are unaffected by the long-range screen- 
ing, at least in pure systems, while the constant term 
T^f (oji) is screened in the same way as the monopole 
charge q cc (k, k + q) » 1 [H H| H3] • 

The Raman spectra, associated with imaginary part 
of Eq. |2T|) . comprise the incoherent electron- hole con- 
tributions characterized by the cut-off frequency of the 
order of qvp and, for the Ai g channel, by the plasmon 
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FIG. 3: Two typical quadratic (a) and bi-linear (b) direct Ra- 
man scattering processes in the conduction band proportional 
to (-Hi) 2 . The self-energy parts on the diagrams treated as 
constant are encircled [41| . The crosses represent the quasi- 
elastic scattering H[. 
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contribution related to the screening of 7^ {iOi). These 
spectra are directly related to the dynamical structure 
factor S(q,oj) — — Im{xi : i(q, u)}. The intensity of both 
the collective and incoherent electron-hole contributions 
to — Imlx^^q, lj, uji)} is proportional to small q 2 . These 
types of signals have never been detected in the high-T c 
cuprates [U, [3(| , in contrast to the semiconducting sys- 
tems, such as GaAs (qvF ~ 50 cm -1 ) 53]. In the high-T c 
cuprates, the measured Raman spectra are roughly pro- 
portional to the optical conductivity, with the intensity 
proportional to the channel-dependent relaxation rates. 
This leads us to study the scattering of the quasi-particles 
on the disorder. 



IV. RAMAN CORRELATION FUNCTIONS IN 
SYSTEMS WITH DISORDER 

A. Incoherent scattering 



FIG. 4: The direct (1 -» 2 -» 3) and indirect (forward, 1 -> 
2 — » 4 — *• 5, or backward, 1 — *• 2 — > 4 — > 6) bi-linear Raman 
scattering processes in the conduction band. The solid lines 
represent the three effective fermionic bands (the indices c, P 
and N) for the typical values of the model parameters A 1 ^ = 
0.66 eV and tf d = 0.73 eV [UEl- The energies are measured 
with respect to the energy of the 2p a oxygen orbitals, so that 
the dispersionless nonbonding band is placed at E p = 0. The 
dashed lines are the photon dispersions, and the dot-dashed 
line is the Fermi energy fi = —1.793 eV corresponding to the 
hole doping S = 0.1. 

FIG. 5: A few direct contributions to the Raman correlation 
functions in powers of (Hi) 2 , according to Eqs. JCjjJl, JCSJ 
and (|C4|) . The full rectangle is the effective Raman vertex of 
Fig. 1(b). 



This section deals with the contributions of the inco- 
herent quasi-elastic scattering to the Raman correlation 
functions x„ tV (q,u),u>i), including the Coulomb screening 
effects. The discussion starts from the low order scatter- 
ing on the disorder, continues by the summations to high 
orders and adds the Coulomb screening at the end. In 
this discussion it is convenient to distinguish between the 
direct and indirect processes, as further explained below. 



1. Direct processes 

As illustrated in Fig. 3, for all correlation functions 
considered in this article (charge-charge, current-current 
and Raman correlation functions), the probability for 
the direct electron-hole pair creation is proportional to 
/ c (k) — / c (k + q) and associated with the resonance con- 
dition Hu> w E c (k) — E c (k + q). The corresponding scat- 
tering paths 1 — > 3 and 1^2-^3 are shown in Fig. 4. 
The direct scattering on the disorder can be roughly in- 
corporated in the correlation functions in the standard 
phenomenological way 52]. Alternatively, one can ap- 
ply the gauge-invariant treatment to sum the direct pro- 
cesses shown in Fig. 5 in powers of (H[) 2 . The gauge 



invariance conserves the number of charge carriers in the 
scattering processes [46j. As shown in Appendix C, for 
lu > qv-p, the latter approach gives the unscreened, direct 
charge-charge correlation function (intra- and interband 
contributions) of the form 

a'LI/ker 

x / £ (k)-,Mk + ) 

hu, + mv L a V (k, W ) + e ll , (k, k) - El ' L '^ k+) ' 

(24) 

where q = £) a ' q a 'e a >, n LL = 1, n Lk = 2, r^' (k, w) = 
Im{E^ L '(k,w)} and J B ii ,(k,k + )= J B i (k)- J B i ,(k + ). 

Eq. (I24| can be easily generalized to other cor- 
relation functions. For the quasi-elastic scattering 
Im{E cc (k, u>)} w r^' d (here, the index i = 1, a, and v for 
the charge, current, and Raman vertices, respectively). 
In the dynamical limit, we thus obtain the universal ex- 
pression for the unscreened, direct intraband correlation 
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k+q k+q 

« <:> * <3r--^o + - 

k k 

FIG. 6: The Coulomb screening of the direct (a) and indi- 
rect (b) processes in the Raman response functions in pres- 
ence of the quasi-elastic scattering. The dotted box includes 
the electron-hole self-energy contributions associated with the 
quasi-elastic scattering processes. 
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FIG. 7: Three typical indirect Raman scattering processes 
proportional to (H[) 2 . The first two include the incoherent 
scattering of conduction electrons while the third shows the 
incoherent scattering in the empty band(s) (the notation is 
the same as in Fig. 4). The effective vertices are encircled 
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direct process indirect process 



—^d T2~ n i,jW- ( 25 ) 



Here l (/i) is the effective density of states at the Fermi 
energy given by 



kcr 



and 

i 2 

q) = 1 in Eq 



i (/i) are obtained by re- 

2 



while < jQ (a*), ra£ 

26J) with (jSf(k))\ 
|7^ c (k, Ui)| , and 7^ c (k, o;i)q cc (k-|-q, k), respectively. Fi- 
nally ja C (k) = ?iJ^ c (k)/(eatpj) is the dimensionless cur- 
rent vertex, Eq. (|A5|) . and is the unit cell volume. 
For the electromagnetic fields (i = a, v) the wave vec- 
tor q = J2 a ' Qa'Sa' is perpendicular to the polarization 
of the fields; i.e. q a i — q z for the symmetrized Raman 
vertices in Eq. (|20|) . 

The RPA series for the screened direct contribution 
to the Raman correlation functions is illustrated in 
Fig. 6(a), and is given by inserting the expression (|25() 
into Eq. ([21"]) . As can be easily seen, the intensity of both 
the plasmon and electron-hole incoherent contributions 
to — Im{x^„(q, 0J,Vi)} remains proportional to small q 2 . 
Fig. 6(b) represents the quadrupolar analog of the well- 
known Hopfield series [45|. It will be argued below that 
the latter is not important for the Raman scattering on 
the high-T c superconductors. 



2. Indirect processes 

Omitting again the Coulomb screening to begin with, 
the disorder-assisted, indirect electron-hole contribution 
is associated to / c (k) — / c (k'), with uncorrelated k and 
k' [see the 1 — > 4 — > 5 processes shown in Fig. 4 
and the related diagrams in Fig. 7(a), as well as the 
1 — > 2 — » 4 —> 5 processes represented by the diagram 
in Fig. 7(b)]. These types of processes become impor- 
tant when fiLO 3> \E c (k) — .E c (k-r-q)|, with the resonance 





FIG. 8: (a) The direct and indirect high-frequency contribu- 
tions [proportional to (.Hi) 2 ] to the Raman correlation func- 
tions [Fig. 3 and Fig. 7(a,b)]. (b) A few indirect leading terms 
in powers of (H[) 2 . 



at hut w E c (k) — E c (k'). This is a typical situation en- 
countered in the absorption and/or emission of photons 
by conduction electrons, i.e. in the intraband optical- 
conductivity and Raman experiments on metals. On 
the other hand, the indirect Raman scattering processes 
1 — * 2 — > 7 — > 6, shown in Fig. 7(c), are directly related 
to the indirect interband optical conductivity [52]. For 
the time-dependent H[ they are essential for the Raman 
analysis of the insulating and semiconducting systems [H . 
In the present case, H[ includes only the quasi-elastic 
scattering and therefore the diagram in Fig. 7(c) has the 
resonant behavior similar to the diagram in Fig. 7(b). 
Thus the processes in Fig. 7(c) can be included in the 
effective Raman vertex (fl"9)) and will not be discussed 
hereafter. 

The direct and indirect scattering processes, shown in 
Fig. 8(a), are large in the high-frequency limit [oc (H[) 2 ]. 
The first qualitatively important corrections to the indi- 
rect high-frequency term come from the second and third 
term in Fig. 8(b) which are proportional to (H[) /u>, i.e. 
they are singular in the zero- frequency limit. The consis- 
tent treatment of the indirect Raman scattering processes 
requires thus the summation to infinity of the most singu- 
lar terms in powers of {H^) 2 /ui. This requires summing 
the singular contributions to all orders in (iJ{) 2 /w in or- 
der to obtain the description which is correct in both the 
high- and low- frequency limits. 

As explained in Ref. [41] in the example of optical 
conductivity, the gauge-invariant treatment of the single- 
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(a) 





(d) 



*1 



FIG. 9: (a) The effective Raman vertex (open rectangle) 
in the indirect processes, [7" (k, uii) — 7£ c (k', Wi)]V r 1 cc (k — 
k')/(ftoj). (b) The expansion of the indirect contribution to 
the Raman correlation functions in powers of (H[) 2 /u, with 
the leading term explicitly shown in (c). The shaded box is 
the electron-hole propagator which is obtained by the self- 
consistent solution of the equation shown in (d) |4l|. The 
diamond is the electron-hole self-energy containing both the 
single-particle self-energy and vertex corrections. 



particle self-energy and vertex corrections in the indirect 
processes gives rise to effective vertices in which there 
is a complete cancellation of the scattering processes as- 
sociated with the constant terms in the bare vertices. 
In the case of optical conductivity, this means that the 
indirect processes in the charge-charge correlation func- 
tions are absent altogether because the effective vertex 
[g cc (k,k) - ^(k'.kOJV^k - k')/(M vanishes due to 
the fact that q cc (k, k) ks 1. In the Raman case, the ef- 
fective vertices ^(k,^) - 7^ (k', w i )]F 1 cc (k - k')/(M 
in Xvvi^i^i) are gi ven by the sum of two terms shown 
in Fig. 9(a), setting 7= c (k,Wi) = 7^ c (u; 4 ) + 7^ c (k, w,), and 
reduce to [% c {k,uj,) - % c {k', Ui)]Vf°(k - k')/(M- The 
contribution to Xuui^^i) °f the constant terms 7^(0;,), 
present only in the v = A\ g channel, thus vanishes, in 
analogy with the case of optical conductivity. In this 
way, xj7 I/ (cJ,o>i) = Xuui^^i) with the hat again indi- 
cating that only the dispersive terms 7^ c (k,Wj) in the 
Raman vertices contribute to Xvvi^i^i)- 



Turning now to the evaluation of x5, d „(ci;, Ui), we first 
note that the leading high-frequency contribution to 
x£\,(o>, LOi) consists of the two self-energy and two vertex- 
correction terms shown in Fig. 9(c). The summation of 
the most singular diagrams in powers of (H[) 2 /oj can 
be performed by using the self-consistent form of the 
electron- hole propagator [4lj . as illustrated in Fig. 9(d). 



(For more details see Ref. |41(.) This approach gives 



,1 Z-^i 



df c (k) (|n cc (k-kQ| 2 ) 

kk'cr Cy ' 

x 7 ^(k,^)(7^(k,^)-7r(k',^)) 
xi[P cc (k,k',c)+P cc (k',k, W )], (27) 

where (l/fr)£>g c (k, k', to) is a useful abbreviation for 

1 



hw + E c (k) - E c (k') + ihr) ' 



Here (...) denotes averaging over the momentum trans- 
fer by the disorder. ££ c (k, to) is the channel-dependent 
electron-hole self-energy, 



^-(k, W ) 



■]T<|vr(q'-k)| 2 )(i 

q' 
1 



h 



tO 



7j c (k,w,) 
Dg c (q',k,c)]. (28) 



The result of the summation of diagrams in powers of 
{H[) 2 /uj is thus 



71 L ' 1 



ka 



1 dE c (k) lu 



-S= c (k,w) 



-E-(k,w) 



7r(k,^) 



2 0/ c (k) -S° c (k,o;) 
0£7 c (k) w + E°°(k,uO 



(29) 



This result has the correct limit for small Hu) in compar- 
ison with typical damping energies. 

It is important to realize here that the expression (j2"9")l 
is obtained under an assumption that is valid for the 
quasi-elastic scattering processes, namely that the real 
part of S£ c (k, to) is negligibly small. In this case, we 
can write S£ c (k,w) « ilm{££ c (k, w)} « ir^ id . This can 
be easily generalized to weakly inelastic incoherent scat- 
tering by introducing T^ ld (w). On the other hand, the 
introduction of -Haf, Eq. J5J), leads to large coherent ef- 
fects in Re{££ c (k, uj)}. This requires the re-examination 
of the single-particle Hamiltonian Hq + Haf, with those 
coherence effects related to Haf incorporated also in new 
effective vertices, and not only in Re{££ c (k, uj)}. The de- 
scription of this procedure is postponed to Sec. VIC. 

The generalization to other correlation functions gives 
the universal expression 



-ir 



c,id 



1 



ir ^ id v 



(30) 



for £? G (k, lu) w ir^' ld , i — l,a,v. Here nj^/z,^) is the 
effective channel-dependent density of states at the Fermi 
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energy of the form 



i£l^(k, Wi ) 

ktr 



2 0/ c (k) 

as c (k)' 



(31) 



and n' 1 d 1 (/i) and n 1 ^ a (fi) are obtained by replacing 
|7= c (k,Wi)| 2 in Eq. ijSTJ with (<f c (k,k)) 2 = and 
(j£ c (k)) 2 = (j^ c (k)) 2 , respectively. Also, we define the 
related effective densities ri^QijUi) and n^(//,w,) using 
the total vertices and the constant part of vertices in- 
stead of |7^ c (k, a;;) I 2 in Eq. §T$. Evidently rf d = and 



Xi d i(w) = 0. Also, 



1,1 



(/'■) 



and r^' d 



— pc,id 



Both those results are required by the continuity equa- 
tion and the gauge invariance of the intraband optical 
conductivity [4l|. 

Let us finally mention the Coulomb screening problem. 
The effects of the Coulomb forces on the indirect pro- 
cesses are described by the Hopfield series of diagrams 
shown in Fig. 6(b), which is an analog of the Hopfield 
series studied in the context of the optical conductivity 
[4ll [45|. This series is free of the q~ 2 singularity and, 
for a sufficiently large relaxation rate r°' ld (with the crit- 
ical relaxation rate defined precisely in the follow- 
ing subsection), does not affect the spectra in a critical 
manner. Therefore, these corrections [starting with the 
second term in Fig. 6(b)] are neglected in the present 
analysis, i.e. we take « x^ d (tJ,Wi). 



B. Direct vs indirect contributions 

When the direct and indirect processes are combined, 
we obtain the total Raman correlation function in the 
form 



X^T(q ) w,w i ) « xS,*(q,w,Wi) (32) 
where 



+xl!.i(q:W,w l )^— rXi„(q,o;,Wi), (33) 

using again the separation of vertices 7^ c (k, u>i) — 
7£ c (o>i) + 7^ c (k, LUi) and the corresponding separation 
of x d 3 (q, uj,uii). There is a well-defined exclusion rule 
here. The constant terms in the vertices participate in 
the direct processes and are screened by the long-range 
Coulomb forces. On the contrary, only the dispersive 
terms participate in the indirect processes. They are in- 
dependent of the wave vector q and are thus nearly unaf- 
fected by the long-range screening. The intensity of the 
former processes is proportional to small q 2 , except in the 
static metallic limit, and the intensity of the latter pro- 
cess is proportional to the channel-dependent relaxation 
rates T^ id . 

To find out which of these two processes dominate 
the correlation functions of the high-T c cuprates, we 



roughly by hT°A 



a ltpd> with aq 



now compare the imaginary parts of the expressions (|25|) 
and O. For a*) « a£* (/*,<*) and T^ d « 

r^' ld , we obtain the condition Huj ps aqtf^. Further- 
more, — Im{Xv V {w,uii)} is characterized by a maximum 
at u> = rf; ld , and the critical damping energy is given 

« 10~ 3 typi- 
1 eV, the re- 

d /(27rc) w 10 cm -1 . For the usual experi- 
mental geometry in the high-T c cuprates [q a i = q z and 
nt tV {(J>,0Ji) w (ti/^) 2 n^(/j, LUi)}, on the other hand, 

the critical relaxation rate is / (2ttc) aqt±, i.e. well 
below 10 cm -1 . Based on this estimates, for frequen- 
cies of the outmost experimental interest, uj/(2ttc) > 50 
cm , the direct processes can be omitted and we con- 
tinue the analysis with the approximate expression 



cally. For the 3D systems and 
suit is r^Q 



xJS al (q,w,w i )«^(w,Wi) 



-id 



(34) 



The measured Raman spectra — Im{x^ al (w, u;,)} are 
thus proportional to — Im{xj7,,(w, Wj)} of Eq. (j3"0)) for ar- 
bitrary to. 

For comparison with experimental and previous theo- 
retical results, it is useful to rewrite the effective densities 
n^„(/i,Wi) in terms of the related densities n^ d (/i,cJi), 
which involve the total Raman vertices. For this purpose, 
we notice that the constant terms 7^ c (wi), defined by 
n d 1 (/i,w;) = 7^ c (wi)nf can be formally expressed 
in terms of the effective density of states n| d (/x, Wj) in 
the following way 



(35) 



This finally leads to 



«i d iM 



hl(p,LUi) = nl((j,,LJi), v = Big, B 2g (36) 
[using the abbreviation h]^(^,,u>i) = r?„ v {}i, u>i)]. 



C. Comparison with the usual field-theory 
approach 

For the sake of comparison with the common field- 
theory approaches (FTA) it is appropriate to notice that 
X d „(q, uj, uii) of Eq. ([33]) can be rewritten as 



x d ,i(q,^)xf ,(q,w) 



xi,i(q,w) 
x d i(q^)x d ^(q^) 

Xi,i(q,w)e(q,w) 



(37) 



in the simplified notation [u>i is omitted and it is noted 
that Xii(q>w) = Xi,i(qi w )]- The relation (071) is also 
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the starting point of the FTA analyses of the electronic 
Raman scattering [29|, [H, HH, |3{|, and is the source of 
controversies regarding the role of the long-range screen- 
ing in the Raman scattering. 

Most of the FTA Raman analyses [H, H HI] 
use the standard approximation for the transverse cor- 
relation functions [2y, |4f| to study the Raman spectra 
in the B\ g and B 2g channels. In this case, x„,„(q, w) 
equals x^ ^(q, w), with Xv,v{<\, given by the diagram 
of Fig. 2(a) for xtviO-i 00 ) m which the momentum re- 
laxation is replaced by the energy relaxation. Equiva- 
lently, this can be formulated by redefining the single- 
electron Green functions with respect to the Green func- 
tions used in the charge-charge correlation functions 
[26L [4f|. For the scattering on the disorder, this leads 
roughly to XvA<h") = X™M = X'^M " xl d „(0), 
with X^vvi 10 ) given by Eq. (|30p . The same approxima- 
tion was extended to the A\ g channel of the high-T c 
cuprates in Ref. (l6j . This is a reasonable approxima- 
tion for the nearly half-filled conduction band with the 
Raman vertices treated explicitly, because the resulting 
ratio Xa.! ,i( t l) a; )/Xi lfa^) turns out to be negligibly 
small, as shown below in Sec. VB. 

On the other hand, the usual approximate descrip- 
tion of the Raman vertices used in the FTA approaches 
generates xl lsi i(q, w)/x d i(q, u) comparable to unity. 
This induces a quite large constant term in the A\ g Ra- 
man vertex, and, consequently, activates the long-range 
forces, as does our approach for a partially filled con- 
duction band. The FTA approaches combine further the 
Coulomb screening in the expression (|37p with the afore- 
mentioned approximation for the transverse correlation 
functions Xu,v{<l,u). The Coulomb term in Eq. (|37|) is 
first removed on taking [29|, [33| Xi,i(q,w) ~ (vFq) 2 /uj 2 ^ 
i.e. the static screening on the ideal lattice in e(q, w). 
Next, the momentum relaxation in j^ .(q,a;) is replaced 



Xi,j l w J 



x^H- 



in the braces of Eq. (|37j) by the energy relaxation. Again, 
this amounts roughly to the replacement of xfj(q, w) by 

one obtains the 
[H [H for the 
screened Raman correlation function in all three channels 
Xiv/(q,w) « x™M; where the Xij(<L,v) are replaced 
by 



-xj d j(0). In this wa- 
common field-theory expression [H, [2! 



FTA 



( W ) 



(38) 



in the braces of Eq. ([57)1 . At frequencies u> -C w p i, the 
form of the resulting Im{x^ A (w)} is thus quite similar 
to the imaginary part of our expression (|34|) [combined 
with ([30|). (|35j) and ([36])]. The background of this result is 
that the large Coulomb term introduced by FTA for any 
band filling is removed therein by the static screening. 

However, instead of removing the last term in Eq. ([57]) 
by the use of the static screening on the ideal lattice 
[23 . [331 1 . our approach determines explicitly the role of 
the long-range forces in the presence of the disorder for 
the typical Raman regime u> > i>F9a, with the Raman 



vertices treated explicitly. It turns out that the last term 
(negligible for the half-filling) is removed from the Ai g re- 
sponse for the partially filled band by the dynamic, rather 
than by the static screening of the long-range Coulomb 
forces involved in the direct processes. This screening 
is characterized by Xij(<la,u) oc q^/ux 2 in all suscepti- 
bilities appearing in the last term of Eq. (|57|) . In addi- 
tion, our approach shows immediately that for ui > r^'j d 
Eq. (|37[) is valid in the impurity-free form, i.e. that the 
plasmon peak does not appear in the Raman response 
due to the xtjila,^) oc q 2 /uj 2 behavior. In contrast to 
that, the FTA does not give a clear recipe how to extend 
its treatment of the last term in Eq. (J37J to the frequen- 
cies u ss Lu p \. It is noteworthy that if e(q, ui) were to 
be replaced here by the usual plasma expression for the 
impurity- free lattice but the behavior of other xfj (q> w)'s 
in this term was kept constant in the small q limit, us- 
ing the expression (|3"8")) . the observation of the plasmon 
would be predicted in the Raman scattering, with a mag- 
nitude comparable to that of the single-particle term in 
xY,T, (uj). This behavior, common in some semiconduc- 
tors [531 ], does not occur in the high-T c cuprates. 

In summary, the Coulomb screening, instead of be- 
ing all-important in the Raman response of the high-T c 
cuprates is not important at any uj. Eqs. ([55]) and (|30[) . 
although widely used, are thus derived here for the first 
time in a consistent manner for u) < r^'' d and extended 
to frequencies around the intraband plasmon frequency. 



INTRABAND RAMAN SPECTRAL 
FUNCTIONS 



In order to illustrate the importance of the enhance- 
ment of the electronic Raman spectra by the interband 
resonance, we shall consider now the bare correlation 
functions X^i/fa^i), v = Ai g , Bi g , B 2g , in the Drudc 
regime of the Haf = case, using (i) the static-Raman- 
vertex approximation, 7^ c (k, uu- u u.) « 7,f (k) , usual in 
most of the current literature [H, [29|, H(J El, [H, [H, 
[H, HH, [36l, and (ii) the elastic-Raman-vertex approxi- 
mation [39l . [55( 1 . Also, the reduced correlation function 
XA ls ,A lg ( w ^) wiU be compared to XA lg ,A lg ( u > u i) to es " 
timate the reduction effects present in Eq. Since, 
in the numerical calculations discussed below, the 3D na- 
ture of the problem appears only in the relaxation rates, 
which are assumed to be independent of the wave vector 
and frequency, we set t± ~ and replace the 3D integra- 
tions in the correlation functions by 2D integrations. 



A. Intraband (Drude) Raman scattering 

With ££ c (k,a;) « ir£ ,ld , the spectral functions related 
to the Drude part of electronic Raman spectra are given 



12 



by 



;rf; id 



(r^ Id ) 



. , . ni d (M,^). (39) 



For r^ :ld w r c,ld , the three Raman channels are still 
distinguished by the effective Raman density of states 
n„(n, u>i) (whatever is t pp ). 

Furthermore, the comparison with the intraband opti- 
cal conductivity 



pc.id 



™tf d \ 2 1 



h J v 



(40) 



with {eat° p s Jh) 2 n lA {^)/v a = e 2 n°*/m, where is the 
effective number of conduction electrons per unit cell 
(discussed in more detail in Sec. VIA), gives an analog 
of the well-known relation valid in simple Drude metals, 



Im {xU(w,Wi)} oc wReKJw)}- 



(41) 



[Notice that h^(p) = n ld (^i), because the constant term 
in the current vertex is equal to zero, i.e. j a (— k) = 
—j a (k).] Here it applies to the CuC>2 plane (a = x, y 
and v = Ai g , Bi g , Big)- This relation has been verified in 
the measured spectra of the overdoped high-T c cuprates 
0,0, Hi, HQ, where the relaxation rates T^ id and T c / d 
have been replaced by T(ui) « T(0) + \ui. 

The SRVA version of these expressions, which sets 
u>i = in Eq. (|39|). was first derived by Zawadowski and 
Cardona [28| and then extended to the case of strong 
quasi-particle damping (A ^ 0) in Refs. [H, [3(| l34j . 
For the overdoped compounds, the A ^ single com- 
ponent intraband term is the only contribution rele- 
vant to the experimental spectra. On the other hand, 
in the underdoped regime, the complete model includes 
both the Drude contribution (l39l and the contributions 



of the lo w-ly ing excitations across the AF (pseudo)gap 
0, 0, [H, [Ef , discussed further in Sec. VI C. 



B. Static-Raman vertex approximation 

We present now the bare spectra —lm{x 1 ^ iy ((jJ, uji)} rel- 
evant for the Drude regime. First we discuss the valid- 
ity of the SRVA [uJi = in Eq. f3T])], As mentioned 
above, the three-band model used in the present calcu- 
lation includes the site energy splitting = E p — 

and the first-neighbor bond-energy t^, but neglects the 
second- neighbor bond-energy t pp [i~4l.Tl5l. |43|. which re- 
stricts our physical discussion to the La2Cu04 family, 
where t pp does not seem to play an all important role. 

The effective density of states n a ([i) was evaluated pre- 
viously [40] for the parameters required to give a rea- 
sonable agreement with the measured spectral weight 
of the visible conductivity in the La2Cu04-based com- 
pounds. Using the same parameters, n„(/i) is calculated 




-2.5 -2^-1.5 -1 -0.5 0.5 
Fermi energy (eV) 

FIG. 10: The dependence of the effective density of states 
on the Fermi energy fi, for Af d = 0.66 eV, tf d = 0.73 eV 
and tpp — 0. The label v — 1 denotes the ordinary density 
of states (divided by a factor of 5), and v = A\ g , Bi g , and 
£>2 9 correspond to three Raman polarizations. For clarity the 
E>2g density of states is multiplied by 10. The hole picture is 
used, i.e. the upper band boundary corresponds to the hole 
doping 5 = 1 (measured with respect to half-filling). The 
doping range < 5 < 0.3, relevant to the hole doped high-T c 
cuprates, is indicated by two arrows. 



now in the SRVA. Fig. 10 shows this effective density 
of states, representing an appropriate measure for both 
the maxima in the Drude part of the Raman spectra, 
Eq. (j3"9"|) , and the corresponding spectral weights. The 
most striking result is that in the doping range of in- 
terest, < S < 0.3, the ratio nB lg {y)/nA lg {lj) is large 
[typically nB lg {^)/nA lg (^) ~ 50]. This enhancement is 
related to the fact that, for t pp = 0, the factor [7^ (k)] 2 
becomes negligible in comparison with [y B c (k)] 2 for the 
Fermi energy close to the van Hove energy. This predic- 
tion of SRVA is however physically unacceptable, since 
the measured n Bls (/i)/nA ls (/i) ~ 1 [ESQ- 

Using the definition of the constant terms 7~Jf(a;i), 
n d uji) = 7^ c (w;)n d we can write 



nl d ^)=nl d (p) 



(M) 



<i(m) 



(42) 



in the simplified notation {u>i is omitted). This expres- 
sion [and its approximate version (|36p as well] reveals 
the existence of two qualitatively different regimes: (i) 
For the nearly half-filled conduction band, i.e. for the 
Fermi energy close to the van Hove energy, the second 
term is negligible [n d x (^) crosses zero at 5 « 0.3, and 

n ii(/ i ) ls singular for /i « £ v h]- For /1 rj £ v h, the con- 
stant term in the A\ g Raman vertex is negligibly small, 
(ii) On the contrary, for the doping well away from half- 
filling, the dispersive terms in the vertices are negligible, 
leading to the strong reduction effects in Eq. (|42|) with 

<>)<<<>)■ 
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FIG. 11: Inset: The resonant enhancement of the Ai g den- 
sity of states riA lg for 8 = 0.1 (p = —1.793 eV in Fig. 10), and 
ftr intc r = 0.1 eV (solid and dotted line) and 0.15 eV (dashed 
line). The dotted (solid, dashed) line represents the contribu- 
tions of the real (real and imaginary) part(s) in 7^ c (k,^i) to 
riA lg ■ Main figure: The total bare effective density of states 
for all three Raman channels for hT intcI = 0.1 eV. The B 2g 
spectrum is again multiplied by 10. 

To simplify the discussion of the resonant effects and 
the effects of the AF correlations, in the rest of the ar- 
ticle we consider the effective density of states n|f (//) « 
njf (/x). For < S < 0.3, the corrections are of the order 
of few percent, i.e. they are comparable to the effects 
of the orthorhombic distortion on xj,°£ al (q, u), LOi) which 
have been already neglected here. 

C. Elastic- Raman vertex approximation 

We calculate therefore the effective density of states 
n v {jx,u)i) in the ERVA, i.e. retaining tOi in Eq. (|3Tj) . 
for the hole doping S = 0.1 and the damping energies 

ftpnter = q j and q e y^ j n p ig n we show thg results 

for the model parameters used above (A°* = 0.66 eV 

and = 0.73 eV). For fkui w 0, the large nB lg {p,^i) 
intraband term, associated with van Hove singulari- 
ties, is large with respect to the interband n,A lg {[i, oji) 
term. For hiOi around i?jv(k) — fj, « 1.8 eV (N for 
the non-bonding band), the resonant (interband) con- 
tribution to nA ls (/i,Wi) is nearly equal to the sum of 
the static (intraband) and resonant (interband) terms in 
riBx In the maximum, the comparable interband 

contributions dominate. This energy range corresponds 
to E N (k) - [i < huji < (E N (k) - £? c (k)) max , because 
for i pp = 0, the optical excitations between the conduc- 
tion and antibonding bands are negligible [43j]. For t pp 
large enough, the latter excitations become important 
as well, and resonant effects are extended to the energy 
region E N (k) - /i < fkJi < (E P (k) - B c (k)) max (i.e. 
between 1.7 and 4 eV). Due to the resonant enhance- 



ment of the Raman scattering processes, we find the ratio 
n B lg (^^i)/nA lg (fJ.,tOi) consistent with the experimental 
observation. Notice, however, the reduction of the reso- 
nant effect with increasing damping energy fir intcr (inset 
of the figure). 

The spectral weight of the Big channel relative to two 
other channels turns out to be one order of magnitude 
smaller than the one usually found in experiments. This 
reflects the fact that various processes described by other 
parameters of the three-band model, and in particular by 
the direct oxygen-oxygen hopping i pp , are absent here. It 
should be noticed that t pp opens an additional channel 
in the electron-photon coupling [see Eq. (|A3|) ] involving 
predominantly the electronic states in the nodal k x = k y 
region of the Fermi surface. As easily seen, this leads 
in the first place to the enhancement of the B2 g Raman 
spectra giving the contributions proportional to t pp in 
7|j c 2 (k,Wi), additional to the contributions of the indi- 
rect oxygen-oxygen hopping processes [oc (t^) 2 ] shown 
in Figs. 10-11. 

We notice finally that, if the contributions of 
Im^^^k, LOi)} to n,y(fi,uJi) are neglected, one obtains 
the resonant structure characterized by two peaks split 
approximately by the energy 2fiF lnter , as represented in 
the inset of Fig. 11 by the dotted line. Similar de- 
pendence of the Raman spectra on the photon frequen- 
cies was already proposed in the multiband study of 
the electron-mediated photon-phonon coupling functions 

It should be noticed that most of the recent Raman 
studies are focussed only on the B\ g and B^g channels. 
These two channels scan the complementary parts of the 
Fermi surface (the vicinity of the van Hove points in B± g 
and the nodal region of the Brillouin zone in Big) and al- 
most all relevant physics is present in the related spectra 
0i Ell Hil ■ Our comparison with the experimental data, 
given in Sec. VI C, will be thus also limited to these two 
channels. 



VI. EFFECTS OF THE AF ORDERING 

In order to make our analysis of the coherence factors 
analytically tractable we shall restrict it here to the situ- 
ations in which the direct oxygen-oxygen hopping t pp is 
not qualitatively important and set it equal to zero. Such 
is the case of La2Cu04 based families for the doping not 
too far from the optimal doping, where the Fermi surface 
is nearly square. 

The effective AF potential A(k) is assumed to be of the 
d x 2_ y 2 symmetry, A(k) = 0.5AAF(cosk • ai — cosk ■ &2). 
This potential dominantly affects the states close to the 
van Hove points, leads to the dimerization of the bands, 
and is accompanied by the low-lying interband processes 
characterized by a threshold energy proportional to the 
magnitude Aaf- Two subbands of the conduction band 
will be denoted by the indices L = C (upper band) and 
L = C_ (lower band). For the half-filled conduction band 



14 



of the t T 



model, Qaf leads to the ideal nesting of 



the Fermi surface and, correspondingly, the relevance of 
this perturbation grows with decreasing hole doping. 

In other cases, the interplay between t pp and A(k) 
is probably responsible for the anomalies regarding e.g . 
the development of both the Fermi surface shape [23| 
and the optical conductivity with doping. Namely, for 
tpp large enough with respect to tf^, even small changes 
in the hole doping could produce dramatic changes in 
the electrodynamic features of the electron system (this 
might be analogous to the situation found in the quasi- 
one-dimensional Bechgaard salts [59] ) . ARPES measure- 
ments in the YP^CaOr-^ and Bi-based cuprates [HI, [HI 
are indicative of such a regime, not discussed here. 
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A. Hall coefficient 

In the three-band model with the magnetic field nor- 
mal to the conduction plane, the room-temperature Hall 
coefficient is of the form i?u ~ l/(ecnn), where nn is the 
effective Hall number given by n H = nf x nf y /nf y . The 
diagonal and off- diag onal effective numbers of charge car- 
riers read as [4(| l6Ctl6l| 

ff _ ml v . cc floa a/cfr) , , 

n %i = ^2y^2 dEc(k) J ^ C ( k )[ 7 «/ C ( k ) J ^ C ( k ) 

k* a 

- 7 ^ c (k)J a cc (k)] (44) 

(a = x or y). The structure of the intraband current 
vertices, J^ c (k) [J^ c (k)], and the static Raman vertices, 
7«/F(k) [7^(k)], for the A AF ^ (A AF = 0) case is 
determined in Appendix B (A). For Aaf ^ (Aaf = 
0), k* (k) refers to the new (old) Brillouin zone. The 
DC conductivity can be scaled by the diagonal effective 
numbers, as well, according to the relations (|4"0")) and (|4"5| . 

The effective numbers (j43|) and (|44|) are extremely sen- 
sitive to the correlation effects. In order to illustrate this 
dependence, the numbers n x s x and nn are calculated with 
and without the potential A(k) of the d x 2_ y 2 symmetry 
and are compared to the experimental observations in 
La2- K Sr 2 ,Cu04 [J 52] showing that (i) the change of the 
sign of nn occurs nearly at 5q w 0.25; (ii) nn °c S in 
the underdoped compounds; and (iii) oc i5 for 5^0. 
The results are given in Fig. 12 for Aaf = and 50 meV. 
The main figure illustrates the well-known fact that for a 
pair of bonding and antibonding bands the critical dop- 
ing do, which separates the electron-like doping region(s) 
from the hole-like one(s), is shifted for finite ip^/Ap^ 
{tp P = 0) from 6 = in the positive (negative) direction 
for the lower (upper) band, breaking in this way a simple 
electron-hole symmetry in each of these two bands. For 
the wide conduction band, characterized by A°^ = 0.66 



FIG. 12: Main frame: Effective numbers n% x and |tih| (rep- 
resenting also the DC conductivity and the inverse Hall coef- 
ficient, scaled by e 2 /(mvoT%' ld ) and ec/vo, respectively) as a 
function of the doping level for Aaf = 0. Inset: The effect of 
the AF correlations on n^, and |«h| for the d x i_ y i symmetry 
perturbation A(k) with Aaf = 50 meV, in the hole-doped 
region. The critical doping So, where 11°* = 0, is labeled by 
arrows, n and p denote, respectively, the region of electron- 
like (riH,w|y < 0) and hole- like (riH,W™ > 0) behavior of the 
charge carriers. 



eV and = 0.73 eV, this results in 5q « 0.27, in agree- 
ment with the observation (i). The measured linear in- 
dependences of ri°j (iii) and nn (ii) can be related to the 
mid-infrared (MIR) gap structure, as seen from the inset 
of Fig. 12. It should also be noticed that, for Aaf not 
too large, the position of 5o is only slightly dependent on 
Aaf- More importantly, due to the doubled number of 
zeros of d 2 Ec(k)/dk a dkp (which appear above and be- 
low the original van Hove energy £ v h), the effective num- 
ber rf xy has two zeros, resulting in an additional critical 
doping within the electron-doped range. In crude terms, 
this restores the electron-hole symmetry of the phase di- 
agram of the high-T c cuprates, which is seen in the Hall 
coefficient measurements [lL [421 . 



B. Optical conductivity 

The dependence of the low-frequency conductivity on 
the symmetry and magnitude of the dimerization poten- 
tial A(k) is analyzed in detail in Refs. 0, EH. For the 
sake of completeness we enumerate here the most impor- 
tant results. The two-component Aaf 7^ intraband 
conductivity reads 

2 eff 

intra/ ,\ r 1 e n aa 10 j- ■ ,MIKf, ,\ 

w m UJ + lTa 

(45) 
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FIG. 13: The optical conductivity (|45[) for the anisotropic-s 
potential A(k) = A A p[0.5 + 0.125(cosk ■ a x - cosk ■ a 2 ) 2 ] 1/2 
with A AP = 45 meV, Af d = 0.66 eV, tf d = 0.73 eV, 8 = 0.1, 
Ci = 0.18 and C2 = 0.4. Main figure: HF^ = 30 meV 
and hr™ m = 50 meV (suitable to T = 200 K spectra in 
the La2Cu04 based compounds). Inset: ftr„ — 15 meV 
and hr™ m = 25 meV (T » 100 K). The data measured in 
La2Cu04.r2 at T = 200 K [56( ] connected by the dotted line 
are given for comparison. 



with the effective number of conduction electrons, n®*, 
and the MIR polarizability, a^ R (w), given by 



-$>£?(k)[l-/c(k)], 
v * — ' 

k*er 

1 1 ^ (M 2 |j£ £ (k)| 2 



(46) 



EhcM 



k*cr 

2^ C c(k)[/c(k)-l] 

(^ + mrMiR) 2 - jB 2 c(k) - 



(47) 



The renormalization factors £1 and (2 in Eq. (T4"5j) serve 
here to model the effects of fluctuations of auxiliary 
bosons on the low- frequency optical excitations [4(| • The 
vertex J a — (k) and the energy difference Eccj^) are 
given in Appendix B. 

Fig. 13 illustrates the typical low-frequency spectra 
measured in La2Cu04.i2, compared to the model pre- 
dictions. In spite of its simplicity the model (|45)) - ((47|) 



with t pp = can explain why the MIR structure in 
La2Cu04.i2 is nearly independent of temperature [56j |. 
Namely, at temperatures below the room temperature, 
the position of the MIR maximum hujyiiK ~ 90 meV is 
well above the relaxation rate HT^ 1R and correspondingly 
Swmir ~ 2Aaf, independent of hT^ m . This situation 
strongly contrasts with those observed in the Bechgaard 
salts |59| or in Bi 2 SrCu06 [13] where small Drude spec- 
tral weights (i.e. v n c ^ a <C 1) reveal the interplay be- 
tween i pp (or £b in the Bechgaard salts) and the energy 
scale 2A A f 113 ■ 
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FIG. 14: The B\ g (a) and £>2 9 (b) electronic Raman spectra 
obtained by ERVA for the d x 2_ y 2 symmetry potential A(k). 



The parameters are A° d 



0.66 eV, tf d 



0.73 eV, S = 0.1, 



hut = 2 eV, KTf m = 50 meV and ftF mtcr = 0.1 eV. The 
curves A (B): A AF = (45) meV and ftr=' id = 30 meV. The 
curves C: A = 45 meV and hT c v ' id = 15 meV [with the Drude 
(dotted line) and MIR (dashed line) contributions indicated 
as well]. The B% g spectrum is multiplied by 10. 



C. Bi g and £>2 9 Raman spectra 

Next, we extend the discussion of the AF effects to 
the electronic Raman spectra. In the hole-doped regime, 
the Drude- like contributions and the low-lying transitions 
through the AF (pseudo)gap are given by Eq. (J2U) and 
by 

-Im{x^ R (u,,u,x)} « 1 J2 |7^(k,^i)| 2 [/c(k) - 1] 



k*<T 



xlm 



, (48) 



respectively. Neglecting the effects of A(k) on the in- 
termediate interband processes and applying the static 
approximation for the low-frequency part of the Raman 
vertex, the elastic Raman vertices in the expressions (|39|) 



and (|48p calculated at t pp = are given by 

7 ^(k,^)«7r(k,^)cos 2 ^ 

+ 7 -(k±Q AF ,^)sin 2 ^, 

7 ^(k,^) « i[ 7 ^ c (k,^) - 7 -(k ± Q AF ,^)] sin^(k). 

(49) 
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if(k) is an auxiliary phase defined in Appendix B. 

Again, for A AF < t^ s d ,A^, the d x 2_ y2 symmetry of 
A(k) causes significant effects in the Raman spectra only 
for relatively small doping (i5 < 0.15) when the Fermi 
energy fi is close to the van Hove singularities. The most 
important qualitative results are illustrated in Fig. 14 for 
S = 0.1 and A AF = 0, 45 meV. 

First of all, we observe in Fig. 14 that the MIR peak 
in the optical conductivity is accompanied by a similar 
peak in the Raman spectra, but only in the B\ g channel. 
As a result, the Raman spectral density increases with 
frequency towards a maximum in the B\ g channel, at 
hio w 2Aaf, in contrast to the B 2g channel, where it 
decreases immediately after the frequency Hcu rj ?T£ ,ld . 
This agrees qualitatively with the Raman experimental 
results. 

Second, the observed [l3| doping-induced weakening 
of the Drude part of the B\ g spectra by one order of 
magnitude with respect to the B 2g spectra below 5 « 
0.15, can be related to the (pseudo)gap features in the 
electron dispersion in the vicinity of the original van Hove 
points. Namely, the B\ g effective density of states at the 
Fermi level (shown in Fig. 10) is strongly suppressed for 
|e V H — m| < Aaf- This contrasts with the B 2g case where 
the spectra come dominantly from the nodal region of 
the Fermi surface, unaffected by A(k). 



VII. CONCLUSION 

The electronic Raman correlation functions have been 
calculated here for the Emery three-band model, using 
the distinction between the direct and indirect scatter- 
ing on the quasi-static disorder. It is shown that there is 
a simple exclusion rule connecting these two scatterings 
and the long-range Coulomb screening. The direct pro- 
cesses concern the constant terms in the vertices. They 
are strongly affected by the long-range screening, and, 
in the dynamic limit, participate in the correlation func- 
tions through the contributions proportional to small q 2 . 
The indirect processes include only the dispersive terms 
in the vertices. They are nearly unaffected by the long- 
range forces, and their contributions to the correlation 
functions are proportional to the channel-dependent re- 
laxation rates. It is shown so that in the high-T c cuprates 
the contributions of the direct processes to the Raman 
correlation functions can be safely neglected. Using the 
elastic approximation for the Raman vertices in two [with 
and without the AF dimerization gap A(k)] analytically 
solvable versions of the t pp = Emery three-band model, 
we show than that the resonant Raman scattering pro- 
cesses remove a large discrepancy between the spectral 
weights of the A\ g and B\ g Raman channels obtained in 
the static approximation for the Raman vertices. The 
resulting spectra agree reasonably well with experimen- 
tal findings. It is also shown that the anomalous MIR 
peak in the optical conductivity, observed in the under- 
doped compounds, is correlated with the corresponding 



structure which appears only in the B\ g Raman channel, 
as well as with the measured linear (^-dependence of the 
Hall number. This relation is explained here in terms 
of the A(k) 7^ AF correlations. On the other hand, 
the A(k) = Emery model used to fit the overall band 
structure, a part of which is seen in the ARPES data 
[HI], leads to different results. Particularly important in 
this respect are Raman selection rules. The small energy 
scales observed in the Raman scattering, just as in the 
ARPES data [23| , are therefore better related to the AF 
correlations within the conduction band than to the low- 
energy interband transitions in the strongly correlated 
A(k) = metallic state. 
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APPENDIX A: THREE-BAND VERTEX 
FUNCTIONS 

The coupling of the vector potential A(r) to the con- 
duction electrons of the Emery three-band model is given 
in the usual way [ZtJ , by replacing the hole creation (and 
annihilation) operators in the bare Hamiltonian Hq by 

ft = jt e ie/(fec)(R n +r,)-A(R n +n) /j^ 

(similar for l na ). Here R„ and r; are, respectively, the 
Bravais lattice vector and the position in the primitive 
cell of the orbital labeled by the index I. The Taylor 
expansion in the vector potential of Hq to the second 
order leads to 

Ho -H « H^= £ a<(k >q )Zt +q<r ^,(A2) 

ZZ'kqcr 

where 

a 

e 2 m ^ d 2 Hk l 'Ck) , , , 

q'a/3 

In the Bloch representation, H ext is given by the expres- 
sion ([5]), with the vertex functions 

(A4) 
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{a, 13 = x, y). 

The number of channels in the electron-photon cou- 
pling is equal to the number of independent bond ener- 
gies; and ip P in the Emery three-band model for the 
in-plane processes. For the t pp = three-band model, 
one obtains the dimensionlcss in-plane current and bare 
Raman vertices (a = x or y) of the form [4p| 



sin k • a Q 



£P( k ) = t f d ^-^ sink - aa , 

■cN li \ j.cff 1. 1. 

Jx ( k ) = t pd — sm-k-a 2 cos-k-ai, 



ir(k) 



;ff 2w k . 1 



pd" 



sin -k • ai cos -k ■ a 2 , (A5) 



and 

7 ^(k;2) = 5 a , 
respectively, with 

T LL' i 



m &f d UkVk . 2 1 



sin -k • a Q , (A6) 

Tflxx tk ^ 



eat 



eff 

J--(k) = l^i^'(k). 



(A7) 



Mk, «k, and tk are the auxiliary functions defined in Ref. 
[3, and m xx = h 2 Af d /(2a 2 (tf d ) 2 ) is the in-plane mass 
scale (| ai| = |a 2 | = a). 



APPENDIX B: VERTEX FUNCTIONS WITH AF 

The AF dimerization of the conduction band E c (k) 
caused by Haf is solved elsewhere [!o|. Apparently, 
Haf can also describe dimerizations other than AF (spin- 
Pierls, charge-density waves) . That is, there is no explicit 
spin-dependence in the dispersions of the bands in this 
Appendix. 

The vertex functions important for the present analysis 
can be shown in terms of the auxiliary phase defined by 



tany>(k) 



2A(k) 



£ c (k)-£ c (k±Q A F) 



(Bl) 



The static Raman vertex and the current vertices relevant 
to both the effective numbers (|43l) - (|44| and the optical 
conductivity (|45p are given, respectively, by 



7L°(k) = 7^(k)cos 



2 <^( k ) 



7S c a (k±QAF)sin 



2 ^( k ) 



m 2|jf g (k)| 2 
E C c(k) ' 



and 

J Q CC (k) 
^ C£ (k) 



J^(k)cos 



2 V(k) 



J^(k ±Qaf) shr- 



ink) -J^ c (k±Q AF )] sin ^(k). 



2 

(B2) 

^(k) 
2 ' 

(B3) 



Here E C c{k) = E c (k) - Eg(k) and 
£ c ,c(k) = i[£ c (k)+£ c (k±Q AF )] 



(B4) 



±\/-[E c (k) - E c (k ± Qaf)] 2 + A 2 (k). 



Similarly, the approximate expressions for the total Ra- 
man vertices are given by the expressions (|49p . 



APPENDIX C: LONGITUDINAL RESPONSE 
THEORY IN MULTIBAND MODELS 

We consider the Hamiltonian ^ with H' 2 = and H ext 
given by Eq. (fTTj) . H[ includes only the quasi-elastic 
scattering processes on the disorder. We introduce the 
retarded electron-hole propag ator D LL '(k,k+,k^,k',t) 
defined by (hereafter q = q a e a ) 

P Li '(k,k+,k' + ,k',i) (CI) 

= -m){[Li^t)L>^ a {t),L'i + ^MLvM]), 

and the related induced density 

Sn LL ' (k, k+ , lj) ee Sn LL ' (k) (C2) 
= E \ VLL ' ( k ' k +, k' + , k', u;) q L ' L (k' + , k')^ cxt (q, lo). 



The equation of motion for V LL (k, k + , k^_, k', t) can be 
set into a form analogous to the Landau equation 

[fiw + E L {k) - E L ,(k + )]5n LL '(k) 

= [f L (k) - f L ,(k + )]q L ' L (k + ,k)V cxt (q,L0) 
-it0m{E^ L '(k,uj)}6n LL '(k), (C3) 

where Sh LL (k) is the contribution to Sn LL (k) which is 
proportional to J% L (k) and 

KS?' (k, W )«-^|Fi(q')| 2 J [P iL ' (k, k + + q', W ) 



+D LL '(k + q',k+, W )](l 



J^ L (k + qQ 
^' L (k) 



(C4) 



is the electron-hole self-energy for the case Vj LL (q / ) 
Vi(q'), and 



X>o (k,k',w) 



fej + Si(k) -£z/(k') + i%' 



(C5) 



In expression (|C3[) the fact that the real part of the 
electron-hole self-energy is negligible for the quasi-elastic 
scattering on disorder is taken into account. 

The total induced density Sn LL (k) consists of the in- 
duced charge and current densities [denoted by 5n,Q l (k) 
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and bn\ L (k) [46|]. satisfying the (intraband) continuity 
equation hu)5n% L (k) + E LL (k,k + )5nf L (k) = 0. The so- 
lution of the Landau equation (|C3j) . together with the 
definition for the total optical conductivity 



(C6) 



LL'ker 

and with the relation 



q LL (k + ,k)V ex \q,UJ 



hJ^ L (k) . , 



£i'L(k + ,k) 



iIC>), (C7) 



[corresponding to Eq. (jT2|) combined with the relation 
q a V cxt (c[,uj) = LE° xt (w)] gives 



/ L (k)-/^(k+) 



J^'(k)| 



fiw + iHTLL'(k,cj) + E LL ,{k, k) 



(C8) 



Here till = 1 in the intraband channel, till = 2 
in the interband channel, T^ L '(k,uj) = Im{E^ L '(k, w)} 
and £J iz /(k,k') = £ L (k) - S L -(k'). The related long- 
wavelength susceptibility and the dielectric function be- 



e 2 Xi,i(q,w) = -V— <J aa {uj), 

*■ — ' UJ 



e(q,w) = l + -^yVtr aa (w), ( C9 ) 



with q = ^ o <j e a . The expressions (|C8|) — (|C9[) are the 
generalization of the well-known single-band Landau re- 
sponse functions [5f3]. Obviously, to obtain Eqs. (f22|) - 
(P?5)) of the main text we have to include the contribu- 
tions beyond the three-band model, as well, by adding 
£oo(q, u>) — 1 to the above expression for e(q, ui). 
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